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Abstract 

We calculate exotic cooling rates of stars due to photo-production 
of light particles (mostly scalars and pseudoscalars) which originate 
from the hidden sector of no scale supergravity theories. Using this 
we can restrict the gravitino mass 1713/2- The range of eliminated val- 
ues of 7713/2 stretches over six orders of magnitude and is given by 

2 x 10 8 < — g — < 6 x 10 13 , irig being the gluino mass. Combining our 

m 3/2 

result with the earlier analysis from colliders ( — — < 2.7 x 10 14 ) we 

m 3 / 2 

conclude that — — < O(10 8 ) except for a narrow window around 10 14 . 
Together with the current experimental limit on nig and cosmological 
constraints on m 3 / 2 , albeit model dependent, our analysis shows that a 
light gravitino is on the verge of being ruled out. 
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In the last twenty years a lot of effort went into establishing a signal of a theory 
which we believe would replace the standard SU (3)c <8> SU {2) L ®U (l)y model at high 
energies. In the order of the increasing scale these theories are likely to be: global 
supersymmetry [|lj, supersymmetric Grand Unified Theory (GUT) |2| and supergrav- 
ity H which, being non-renormalizable, should in principle be only an approximate 
version of something more fundamental, like superstrings. Whereas we hope that 
the effects of global supersymmetry will manifest themselves around 1 TeV scale 
the physical effects of GUT theories which are testable at present in laboratory ex- 
periments are only a few (to mention the standard example of proton decay or the 
existence of a monopole). Up to now searches for supersymmetric particles and pro- 
ton decay have resulted in lower bounds on their masses and an upper bound on 
proton lifetime M. The presence of supergravity theory at Planck scale justifies the 
soft supersymmetry breaking terms in the global theory i.e. effectively reduces the 
number of free parameters which via renormalization group equation can be scaled 
down to present energies. A more direct signal of the usual supergravity is practically 
non-existent. Though, in general, not a serious problem, this can be a bit intriguing as 
supergravity theories are not unique due to the arbitrariness of the Kahler potential. 
It has been, however, realized that a superlight gravitino (G) could be the harbinger 
of physics from the hidden sector of supergravity ||. A superlight gravitino field 
i/jp acts as spin-1/2 Goldstino x with ip^ = i^j2/3m^ 2 d^x- I n momentum space the 
spin summed density matrix corresponding to this longitudinal part of ip^ takes the 
form p^Pv/m\i 2 which effectively, when coupled to a gauge boson V (V — 7, g, Z) 

~ TTl " ~ 

and a gaugino V, enhances the coupling by a factor — — . More precisely, since G is 

m 3 / 2 

the superpartner of graviton it couples to particles with the strength proportional to 
k = \^8ttGn = 4.11 x 10 -19 GeV' 1 (Gn is Newton's gravitational constant). Hence 
the effective net coupling turns out to be Kmy/m 3 / 2 . Later the idea of a superlight 
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gravitino was embedded in the no scale supergravity theories 0. Special attention 
was drawn to the super-Higgs mechanism in these type of theories and it was found 
that the coupling of certain light scalar (<S) and pseudoscalar (V) to gauge bosons is 
also proportional n{rriy /m^/2) 0. These particles {S and V), in the physical spec- 
trum of the theory, are members of a chiral superfield in the hidden sector. Their 
spin-1/2 superpartners are Goldstinos which get eaten up by the gravitino. Needless 
to say that a signal originating from the hidden sector of such a supergravity theory 
would be quite spectacular, more so as the relevant coupling is directly proportional 
to the gravitational constant Gn- Indeed bounds (to be discussed later) on m 3 / 2 have 
been put by examining exotic decay processes and scattering reactions involving S, 
V and G at collider energies P-|TI|. It is worth mentioning here that what makes 
bounds on gravitino mass even more interesting are hints from cosmology that point 
in the direction of either a light gravitino with mass m 3 / 2 < 1 keV or a heavy one 



with m 3 / 2 > 1 TeV [|12],|T3|. In this letter we eliminate the parameter space of m 3 /2 
further by investigating exotic cooling processes of stars like the Sun and Red Gi- 
ants. In particular these additional cooling rates are due to the photon reactions like 
r ye ± — * 7 + S/V (Primakoff process) and 77 — * GG, SS, W . 

The interaction lagrangian relevant for single and two photon processes is [fUjjTH 



4V3Vm 3 /2/ 



2 V 2 

where e is the determinant of the vierbein and T>^ is the covariant derivative expressed 
in terms of the spin connection uj^ Viz. V^ipv = (d^ — i/Au^abCr^ipu- A is the photino 
field and is the mass parameter which enters the full lagrangian in the form of 
the bilinear m^XX. T^ v is the usual electromagnetic stress tensor. In general, one 



can have more than one pair of S and V particles, but for practical purposes it is 
convenient to take just one pair. 

The first two terms in eq.(l) are sufficient to calculate 7e — > 7 + S/V. The 
rest of the interaction terms appears in the calculation of 77 — > GG where we have 
two diagrams with 7 in t- and ^-channel as well as three s-channel diagrams with G, 
S, V in the intermediate state. Note also that since the dimension 5 operators ST -T 
and VT ■ T contain higher order derivatives they are expected to violate tree level 
unitarity at some scale. A number of papers |7||| have been devoted to this subject. 
The violation of tree level unitarity need not be a serious embarrassment of the theory 
since supergravity should not be considered as the last word in the physics of Planck 
scale. We view it rather as an approximation to the latter. Therefore as long as we do 
not stretch the involved energies beyond a critical value (in astrophysical applications 
this is never the case) the results are still trustable. 

At tree level S and V are strictly massless. They can, however, acquire mass 
radiatively, mostly through loops with intermediate gluons. Their mass is estimated 
to be |T| 



m s/v ~ « -pM ^qcd (2) 

Commonly one assumes the equality of gaugino masses at a high unification scale. Via 
renormalization group equation this assumption leads to the following mass relation 



at lower energies [16 



8 G^gTTT, 1 , \ 

rrij ~ - m- g ~ - m- g (3) 

o a s o 

In what follows we will use eq.(3) as a guiding relation between m g and rrij. We 

(Tfl~ \ 
— — can be 
m 3 /2j 

translated into a bound on ms/v an d vice versa. Furthermore, the coupling gs/v-ry 



depends now linearly on m^/p in which case astrophysical limits on nig/v and gs/v-yy 
should not be treated independently, a situation encountered also in axion physics. 

C7Tl~ \ 
— — obtained so far comes from 
m.3/2/ 

e + e~ colliders in the reaction e + e~ — > 7 + nothing realized in the model by e + e~ — > 



jS, 7P ||11|| . A slightly model dependent bound, assuming rriq > m g , can be derived 
in pp reactions studying the events of monojets and dijets with missing px jlOfl . This 



would be the signal of processes like pp — > (75, gP, gg, The outcome of the 

two analyses yields 

(^-) < 1.8 xlO 15 , (^-) < 4.5 xlO 13 (4) 
\m 3/2 J e+e _ \ m V2j pP 

corresponding to ms/v — 0.18 MeV and ms/v — 4.4 keV, respectively (Aq C d = 
200 MeV). Both limits are better than the older Fayet's bound of 4.3 x 10 15 . 

The importance of restricting gravitino mass becomes also clear in light of 
the cosmological arguments given in [|l!|,|13|]. Based on the observed bound on cos- 
mological mass density and a Helium synthesis there are strong indications, quite 
independent of the specifics of an underlying supergravity theory, that 

1 keV m 3/2 <^ 10 TeV (5) 

i.e. cosmology seems to favour either a light or a heavy gravitino. This restriction 
can be evaded by assuming an inflationary scenario with low reheating temperature 



17|. This is, however, incompatible with baryogenesis from GUT theories which 



should occur at higher temperature QT8 1 . A natural way out would be to produce 



baryon asymmetry through anomalous processes at electroweak scale [[nj. Whether 
this is possible remains still an open question, but there are convincing arguments 
which show that the Higgs mass required for such scenario is in contradiction with 



recent LEP limits on Mm ggs [Til. The baryon asymmetry produced at GUT scale 



with B — L consevation gets washed out by the anomalous processes unless e.g. some 
lepton number violating processes are operative ||20|| . An alternative would be a GUT 
theory with no B — L conservation where the baryon asymmetry can, in general, 
survive the electroweak transition. 

Yet another way to avoid restriction (5) can be broadly described as belonging 
to the class of R-parity breaking theories |^T[. In any case, we will envisage a change of 



our understanding of supergravity theories and/or cosmology, should the limit given 
in eq.(5) be violated. 

Before applying astrophysical processes ]22| to restrict the gravitino mass let 



us make some general comments on what is to be expected in such an analysis. From 
the core temperature of the stars like Sun (T ~ 1 keV), Red Giants (T RG ~ 10 keV), 
White Dwarfs (T]y D ~ 2 keV) P3| and a newly born Neutron Star in a Supernova 
explosion |24]] (T S n — 30 MeV) it becomes evident that only light particles (in the 
extreme case with a few MeV mass) can contribute to the cooling processes of stars. 
Plasmon effects will suppress single and two gamma processes (relevant in our case) 
in White Dwarfs and Supernovae. Hence we are left with the Sun and the Red Giants 
where the available energy is at most few keV. Moreover, if the particles S and 
V become too heavy their relatively fast decay will preclude the possibility of their 
causing any exotic cooling. From the width 

2 / \ 4 



1 of m~ 



Y{S/V - 77 ) = — . 2 ( ) m% v ~ rnsp (6) 

96vr \m 3/2 J 1 3456tt \k QCD j 

we can estimate the distance Xfu g ht traveled by the particle in the rest frame of the 

star in meters 

/ keV \ 6 / uj' \ ( ffn s /-p\ 2X 



where u' is the energy of the produced particle. Let us apply eq. (7) to the Primakoff 
process on a target Z = e, p etc. In the case of interest where the initial photon 
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energy c<j 7 is T star we have m z » uj 7 » M s /v- The four momentum transfer q is 

then very low, m S j v y - 2 J < \q \ < 4cj 7 , such that c<j 7 ~ T star ~ a; . If we take 

ms/v ~ 10 _1 T st ar we get for the Red Giants Xfu g ht — 10 16 m. The other extreme 

example of a Supernova with m s /v = 1 MeV (4 keV) gives \fn g ht = 90 (2 x 10 16 ) m. 

The above arguments make it also clear that since we have to ensure that 

the particle takes away energy from the star it cannot be produced almost at 

rest. The analysis becomes then a little bit more involved as we will have to 

start with an assumed value of ms/v, typically ms/v — 10 _1 T S i ar , corresponding to 

f — — j ~ io 12 ( - st< * r r \ , put the limit through exotic cooling on m s /v (alternatively 
\m 3 / 2 J \keV J 

(77l~ \ 
— — ) and finally check if the bound so obtained is better than the starting as- 
m 3/2 J 

sumption. Similar method has been used in p5| for constraining the 77 coupling of a 



light pseudoscalar particle from the physics of pulsars. For consistency reasons there 
one starts with a mass smaller than 10 -10 eV. If the two photon coupling depends 
upon mass, as it is the case for axions in a model dependent way, one has then to 
check, from model to model, if the bound on the particle mass so derived is better 
than IO" 10 eV. 

The first process of interest in our case is the Primakoff reaction 7Z — > 
ZS, ZV with Z = e, p. In the range mz » uj^ » u' one obtains 

^ = a (l Z - ZS) = a( lZ ^ZV) = ^( { 2 h (^-) + 2k2 - l} 

2 6 V m 3/2/ 1 \ms/rj J 



In the approximation we have used the cross section of the Primakoff process is 
independent of the target mass and is the same for spin-0 target (relevant for Helium 
burning stars). Note also that the cross section becomes singular as m s /v —> 0. 

The sum of the cross sections for all possible 77-scatterings in the limit uj 1 > > 
ms/v, ™<3/2 reads @ 



-7 / \ 4 

7 a m~ 



X=SS,VV,SV,GG x •'■ ■ 



The exact calculation of an energy loss (e) per unit mass and unit time involves 
an integration over the initial momenta of the squared matrix element folded with 
the Bose-Einstein or Fermi-Dirac distribution and taking into account the chemical 



composition of the star and possible Pauli blocking pz| . Here we will adopt a sim- 
plified version which is sufficient for our purposes and should give the correct results 
for stars like the Sun and Red Giants. We can write the exotic cooling rates due to 
the single and two photon reactions as 

n e n 7 cr 7 a/ n p n 7 a 7 u)' 

Q Q 
C 77 = 1 U 1 ) 

Q 

where g is the core density of the star and rix is the number density of the particle 
species 7, e, p. The factor 2 appears in eq.(ll) because two particles are emitted. 
With co> 7 ~ to' ~ T and g ~ n v m^ (mjy is the nucleon mass) one obtains 

e 7 * ^a em — f^|V< 2 In [ — I + 21n2 - 1 



Q C-//1 I 

3^ \m 3 / 2/ 
C(3) ™ 2 



urn 

1 V m s/W 



27^-A5^ ^1 { 2,n (^) +21 ■ l2 - 1 } <12) 



„Iffl?! t< KV j.ci^ (13) 

-27^ £> V'»s/J 34992 ^ £>A$ CD 
where £ is Riemann's zeta function. 

(Tfi ~ \ 
— — can be obtained now from eqs.(12) and (13) by impos- 
^3/2/ 

ing e 7 , e 77 < e star . The set of astrophysical parameters needed in the analysis are as 
follows: T ~ 1 keV, g Q ~ 10 2 g cm~ 3 , e & ~ 17.5 erg g~ x sec' 1 and T RG ~ 10 keV, 



Qrg — 10 4 g cm 3 , e RG ~ 10 2 erg g 1 sec 1 . As explained above one assumes a 

starting value, m s /v ~ 10 _1 T star or equivalently f ) ~ 10 12 C^ star \ which is 

\m 3 / 2 J \keV J 

consistent with temperature of the Sun and the decay lifetime of the particle. The 

(Tfl~ \ 
— — < 2.4 x 10 9 and for 
m 3 /2 ; 

— — < 6 x 10 7 . Indeed the analysis is consistent since the 
bound so obtained is orders of magnitude stronger than the starting value. One can 



also perform a slightly better analysis by an iteration process. Translating the bound 

(m~ \ 
— — into the corresponding value of ms/v one can insert the latter back into 
"^3/2 / 

eq.(12) and proceed till a convergence is reached. Note that the convergence is then 
independent whether one starts with ms/v ~ 10 _1 T star . or ms/v — 10 _2 T S f ar . As a 
check one can use the second form of eq.(12) expressed in ms/v- The conclusion is 
then that the following region is excluded 



1.8 x 10 9 i ( i 10 12 , 3.3 x 10 7 i ( i 10 13 



(14) 



These limits justify our analysis as it excludes huge range of the parameter space. On 
the other hand in view of the first (second) bound in eq.(4) it also means that the 

(m~ \ 
— — between 10 13 and 1.8 x 10 15 (10 13 and 4.5 x 10 13 ) remains unrestricted 
m 3 /2/ 

by our analysis. It is worth emphasizing that taking the bound from pp colliders and 
combining it with eq.(14) this window of allowed values becomes quite narrow. 

Let us also make a rough estimate of what is to be expected for a White 
Dwarf (T WD ~ 2 keV, g WD ~ 10 5 g cm~ 3 , e WD ~ 5 x 10~ 3 erg g~ 1 sec~ 1 ). 
Here the plasmon effects play a non-negligible role and can be taken into ac- 
count approximately by multiplying e 7 by the suppression factor e~ m *l T where 
m* is the plasmon mass. As compared to the Sun the gain/loss factor is 



Ti 



i 



WD 



T Q J \e WD J 



( 6 Q ^ e ~m,/T WD 



-1/2 



1.2 x 10 for m^/Ty/D ~ 2. Roughly speaking 
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the bound from White Dwarfs is not expected to be better than the one obtained 
already for Red Giants. 

The two gamma cooling rate is, as one can see from eq.(13), very sensitive 
to the temperature. However, even for stars with high core temperature this cannot 
overcome the other small factors, like k 4 , to make the rate appreciable. For the Sun the 
limit one gets from eq.(13) is m^/p < 3x 10 5 keV which is clearly not a useful bound at 
all. For a Supernova (T S n — 30 MeV, q S n — 10 12 g cm -3 , e S N — 3.6 x 10 19 g cm" 3 ) the 

-1/4 



gain loss factor is 



' T SN \ ( _Qo_\ ( e -2m,/T SN 



10 2 for m*/T WD ~ 5. 



T & / \Qsn) VesTv/ 

This gives ms/v < 1 MeV which is much weaker than the collider bound. Thus the 
Supernova analysis is not likely to provide a better bound than in eq.(14) (this is 
true for the processes we have investigated here, an additional process relevant for 
Supernova could be e + e~ — > 7* — > 7 + S/V). 

(vn~ \ 
— — which excludes 

values over several orders of magnitude. The window around 10 12 — 10 13 which our 
analysis leaves unconstrained can probably be closed in future collider experiments. In 

(7Tl~ \ 
— — < 3.3 x 10 7 could be put. With m„ > 140 GeV 
m 3/2 J 



26| this on the other hand would imply 7713/2 > 0.7 keV which comes close to the 
bound from cosmology (eq.(5)) and hence leaves only a small window for light grav- 
itino. Note that the inherent astrophysical uncertainties might introduce a certain 
margin of fluctuation of the numbers stated above, but the main conclusions remain 
unchanged. 
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